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Abstract 

In a recent publication, we proposed two possible wave functions for the elementary excitations 
of the SU(3) Haldane-Shastry model (HSM), but argued on very general grounds that only one or 
the other can be a valid excitation. Here we provide the explicit details of our calculation proving 
that the wave function describing a coloron excitation which transforms according to representation 
3 under SU(3) rotations if the spins of the original model transform according to representation 
3, is exact. We further provide an explicit construction of the exact color-polarized two-coloron 
eigenstates, and thereby show that colorons are free but that their relative momentum spacings 
are shifted according to fractional statistics with parameter g = 2/3. We evaluate the SU(3) spin 
currents. Finally, we interpret our results within the framework of the asymptotic Bethe Ansatz 
and generalize some of them to the case of SU(n). 

PACS numbers: 75.10.Pq, 02.30.1k, 75.10.Jm, 71.10.Pm 
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I. INTRODUCTION 



In recent years, there has been substantial interest in models in condensed matter physics 
with symmetry groups larger than SU[2), the group underlying the usual spin algebra. In 
particular, transition- metal oxides [l|, Q, S, l^f , where the electron spin is coupled to orbital 



degrees of freedom, have been described theoretically by models with SU(4) symmetry 
{ S 0|- Furthermore, there has been a growing interest in the SU(n) generalization of 



the Hubbard mode 
ultracold atoms 



ll\ motivated by possible experimental realizations in systems of 



mm 



In a previous paper jl4l |. we have used the SU(n) Haldane-Shastry model to make the 
case that the elementary excitations of SU(n) spin chains transform under the representation 
conjugate to the representation of the SU(n) spins on the chain, and only exist in one n th 
of the Brillouin zone. In this article, we will present for one thing the detailed calculations 
underlying our line of argumentation i n |l4| . To this end, we focus on the SU(3) Haldane- 
Shastry model (HSM) Q, Q, H Q, H 3, E , |22j, which serves as a paradigm for not 



only the SU(3) spin chain, but, as we shall see, also illustrates some very general properties 
of SU(n) chains. In particular, we derive the quantum numbers of the elementary, fraction- 
ally quantized excitations, the analogs of the spinon excitations for SU(2), which we call 
colorons for SU(3). As already mentioned, the key result is that these excitations transform 
under the SU(n) representation conjugate to the representation of the original SU(n) spins 
localized at the sites of the chain. In the case of SU(3), if a basis for the original spins 
is spanned by the colors blue, red, and green, a basis for the coloron excitations is given 
by the complementary colors yellow, cyan, and magenta (see Fig. This result is mean- 
ingless for SU(2), as the representations of SU(2) are self-conjugate, but significant in all 
other instances of fractional quantization in SU(n) chains and possible higher dimensional 
liquids, regardless of model specifics. In our analysis, however, we will primarily focus on 
the SU(3) HSM, as this is the simplest model in which this general result can be illustrated 
through exact calculations, and only afterwards generalize our results to the case of SU(n). 
The other significant advancement we report here is the construction of the exact, color- 
polarized two-coloron eigenstates. We obtain those through derivation and solution of a 
Sutherland-type equation which is similar to the case of two spinons in the SU(2) HSM. 
Our explicit calculation implies that the colorons are, just like the spinons, non-interacting 
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or free, but that the spacings in the difference of the individual coloron momenta p n and 
p m with p n > p m are given by p n — p m = ^r(§ + integer) . We will argue that this spacing 
is a direct manifestation of the fractional statistics |23l . l24j of the colorons with statistical 
exclusion parameter g = 2/3. This value is consistent with what we find by naive state 
counting. We then proceed by calculating the SU(3) spin currents, i.e., the eigenvalues of 
the diagonal components of the rapidity operator A 3 and A 8 , for the one- and two-coloron 
states. Thereafter, we interpret our results within the framework of the asymptotic Bethe 
Ansatz. Finally, we generalize some of our conclusions to the case of SU(n). Many details of 
our calculations have been exported, largely in the form of theorems and their proofs, into 
an extensive list of appendices. 



II. SU(3) HALDANE SHASTRY MODEL 

The SU(3) 1/r 2 or Haldane-Shastry model Q, 0| is most conveniently formulated by 
embedding the one-dimensional chain with periodic boundary conditions into the complex 
plane by mapping it onto the unit circle with the SU(3) spins located at complex positions 
r] a = exp(i^a), where N denotes the number of sites and a — 1, . . . , N. The Hamiltonian 
is given by 

#SU(3) = (-jT7j S~] | " 12> (1) 

where J a = \ J2ar c aa^rC aT is the 8-dimensional SU(3) spin vector, A a vector consisting of 
the eight Gell-Mann matrices (see App.[X]), and a and r are SU(3) spin or color indices, which 
take the values blue (b), red (r), or green (g). For all practical purposes, it is convenient to 
express i?su(3) directly in terms of colorflip operators e a J = c^ aa c aT : 

r> 2 N 3 pO-T TO" 1_ 



^su(3) — y. y. 

Ct<f3 <TT y 



EE 



2vr 2 ^ ^ e a J ej tt 2 A^ 2 - 1) 



where the color sum includes terms with a — r. 

The model is fully integrable even for a finite number of sites; the algebra of infinitely 
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FIG. 1: (Color online) Weight diagrams of the SU(3) representations 3 and 3. J 3 and J 8 denote 



the diagonal generators 



25] 



many conserved quantities is generated by the total SU(3) spin 



N 



[H sm ,J] =0, J = J2 J < 



(3) 



a=l 



and rapidity operators 



r tt Aal n a a 1 + VP rabc jb jc 

[if su(3) ,A \ = 0, A = - 2^ - _^ J o'a^, 



(4) 



where a, 6, c = 1, . . . , 8 and / afec are the structure constants of SU(3) defined through 
[A a , A 6 ] = 2f abc \ c . The total SU(3) spin and rapidity operators do not commute mutu- 
ally, 

[j a ,A b ] = f abc A c , (5) 



and generate the infinite dimensional Yangian algebra Y(si3) 18l \26\. 



III. GROUND STATE 



The ground state of i?su(3) for A" = 3M (M integer) is most easily formulated by 
Gutzwiller projection of a filled band (or Slater determinant (SD) state) containing a total 
of N SU(3) particles obeying Fermi statistics (see Fig. |2K) 

i*o> = p c n 4i4 1°) = P G |^D>- (6) 

M<<?F 

The Gutzwiller projector 

N 

P G = []K-2)K-3) (7) 

a=l 

4 




FIG. 2: (Color online) a) Total antisymmetric iV-particle state, b) Typical configuration in l^o)- 



with n a = c^ ah c ah + r c ar + c£, c eliminates configurations with more than one particle 
on any site, and, as the total number of particles equals the total number of sites, thereby 
effectively enforces single occupancy on all sites. As |^so) is an SU(3) singlet by construction 
and Pq commutes with SU(3) rotations, |\l/o) is an SU(3) singlet as well. 

If one interprets the state |0 g ) = Yla=i c ag I ) as a reference state and the colorflip 
operators e bg and e rg as "particle creation operators" , the ground state (jHJ) can be rewritten 



as 
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21] 



(8) 



|*o>= £ *o[z i ;w k ]e b z f...e b / M e%...el* M \O s ), 

{Zi,W k } 

where the sum extends over all possible ways to distribute the positions of the blue particles 
zi, . . . , zm and red particles Wi, . . . ,wm over the iV sites (see Fig. |2t>). The ground state 
wave function is given by 



Mi M 2 Mi M 2 Mi M 2 

i=l k=l i=l k=l 



(9) 



(10) 



i<j k<l 

with Mi = M 2 = M, as derived in App. [Bj The ground state energy is 

18\ N / 

The total momentum, as defined through e tp = ^o[r]iZi, rjiWk\/^o[zi, w k \ with r/i = exp(z^), 
is p = regardless of M. (This is only true for SU(n) with n odd.) 

We will now prove by explicit calculation that Q is an exact eigenstate of the Hamiltonian 
(J2J). Parts of this calculation will be used in Sec. HVBI to prove the exactness of our proposed 
wave function for the coloron excitation. 

To evaluate the action of Hsu (3) 011 © we first replace e| g e| g in (J2J) by (1 — e b — e"){l — 



Hi 



2vr 2 



SU(3) 



N I 



'a °/3 



a °/3 ~ °a °/5 



N 



9-7T 2 1 

zm_ ^ 1 / bb bb , rr rr , bb rr\ 



(11) 



27T 



V - (e hh + e rr ) + -— V - 

and then evaluate each term separately. 

The first term [ek g e| b \I/o] [zi\ Wk), which vanishes unless one of the Zj's is equal to r/ a , yields 
through Taylor expansion 



V P bg gb 



E 



Ml N 

i=l P^i 
Mi N-1 

= EE 



J}1 ^ [. ..,^_i,?7 /3 ,z^ + i,...;wfc] 

-Vp\ 2 



m! dzT 1 2J* 



i=l m=0 

j^(N 2 + %Ml - 6M 1 (N + 1) + 3) 



M a M 2 



Mi 



1=1 fc=l 

Mi M 2 



fc=i 

1 Mi M 2 2 

2 



(12) 
(13) 

(14) 

(15) 



where we have used deg x ,^o[zi] Wk] = JV — 1 and defined A m = — ^2^=1 VaiVa — l) m 2 ■ 
Evaluation of the latter yields A = (N — 1)(JV - 5)/12, A x = -(JV - 3)/2, A 2 = 1, and 
A m = for 2 < m < JV — 1 (see App. EJ). Furthermore, we have used 



x 



+ 



+ 



1, x,y,ze 



(16) 



(z-y)(z-z) (y-x)(y-z) (z-x)(z-y) 
The second term [e^ g eJ| r \I/o] [ z u w k] can be treated in the same way and yields together with 
the first term in (JT3Jl 



Mi M 2 

JV — 3 \ - v z 



Mi M 2 

y^ g ~ J y^ y^ w k 

2 ^ ^ Zi- w k 2 ^ ^ Zi- w k 

i=i fc=i * K «=i fe=i * K 



JV-3 



MiM 2 



and with one part of ()14j) and (|T3j) 



Mi M 2 

EE 



1 

+ 



U7 



~ %')(^ - w k) 2 (z - W k )(Zj - w k ) 



-M 1 (M 1 -1)M 2 , 
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as well as similar expressions for Zi <-> Wk- 
For the third term of (fllj) we obtain 



br rb 



E^^*» 



Ah M 2 Mi / x Af 2 



tr tr (* - V % - * / V m - m J 

Mi M 2 

- <17) 

1=1 fc=l v ' 

Mi M 2 Mi M 2 

wi^fa- z o)i z i -w k ) ° (w k - Zi){w k - Wi) 

Mi M 2 Mi-l , , m _ 2 

+ Z^2_. m \2^f z - Zi )---(z a - Zi ) 1 J 

i=l fc=l m=2 { aj } v ai % > K am % > 

Mi M 2 Mj-l , ]n , ]n _ 2 

Mi M 2 Mi-l M 2 -l / ^„ 



+EE E E 



mini 

i=l k=l m=l n=l 



y- ZjW k (zi - w k ) m+n 2 ^ 

where {a^} ({6z}) is a set of integers between 1 and Mi (M2). The summations run over 
all possible ways to distribute the z a . (w bl ) over the blue (red) coordinates, where Z{ {wk) is 
excluded. The two terms (fT9|) and (}20|) vanish due to 



Theorem 1 Let M > 3, z G C, and zi, . . . , zm G C distinct. Then, 

\M-3 



ii ■ ■ ■ > 

M 



Y, II J ~_ = 0- (22) 

i=l Ylj^i( Z j Z i) 

A proof of Theorem ^ is given in App. [OJ The last term (j21j) can be simplified using a 



theorem due to Ha and Haldane 



2l|: 



Theorem 2 Let {aj} be a set of distinct integers between 1 and Mi, and {bi} a set of 
distinct integers between 1 and M 2 . Then, 

Mi M 2 M1-IM2-I , 1]n 7 ...( 7 ._„..\m+n-2 

\~ <> ' i " k\ - i ~ » /,' j 



m!n! (z ai — ;&) • ■ • (z a — ^Hm, — • • ■ (w b — Wk) 

i=l fc=l m=l n=l {aj}{bi} K 1 7 V ™ 7V 1 7 V m 7 

min(Mi,M2) 

— m). 
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Furthermore, the two terms in line (|18jh together with the remainder of (fPfj) and the cor- 
responding expression from the second term of the Hamiltonian, can be simplified to yield 
M 1 M 2 (M 1 + M 2 - 2)* /2. 

Finally, the diagonal terms of (jllj) cancel the remainder of (|13|) as well as (|17|) and yield 
the additional constant 

1 / \ N 2 -1 /N \" 

-^M 1 {M 1 -l) + M 2 {M 2 -l)}+—^-^--M 1 -M 2 ) . 

When collecting all terms evaluated above and setting Mi = M 2 = N/3 we obtain the 
ground state energy (fTTIJ) . 

IV. ONE-COLORON STATES 
A. Trial wave functions 

We now turn to the heart of our analysis, the elementary and fractionally quantized 
excitations, which we call colorons. In principle, there are two possible, non-equivalent 
constructions for localized excitations starting from @. We may either create a particle 
with color o on a chain with N = 3M + 1 sites before Gutzwiller projection, 

I^HM^Kb 1 )' ( 23 ) 
or annihilate a particle with color o on a chain with N = 3M — 1: 

|^> = ^gc 7(T |^ d +1 >. (24) 

In both cases, or c 7cr creates an inhomogeneity in color and charge before projection. 
The projection once again enforces one particle per site and thereby removes the charge 
inhomogeneity, while it commutes with color and thereby preserves the color inhomogeneity. 
Consequently, the trial states \^ c la ) and describe localized "excitations" of color a 

or complementary color a, respectively, at lattice site r^ 7 . Since i?su(3) is translationally 
invariant, we of course expect neither of them, but only momentum eigenstates constructed 
from them via 

N N 

= N £ e ~^ 7n l*7> = n S^)" ^ ' (25) 

7=1 7=1 

where n is a momentum quantum number, to be eigenstates. 
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27T 2 
iV2 



The important thing to realize now is that only (|23j) or (|24|). but not both, can describe 
a valid excitation. To see this, let us assume a = b for ease in presentation, and note that 
(J23j) is apart from a normalization factor equivalent to 

|^ b ) = P G c 7r c 7g |vl/^ 2 ), (26) 

i.e., creation of a blue particle before projection is tantamount to annihilation of both a red 
and a green particle at site r/ 7 . If momentum eigenstates constructed from via (|2~5|l 

were energy eigenstates, the anti-red (cyan) and anti-green (magenta) coloron excitations in 
(f2Tlj) would individually seek to be momentum eigenstates, which implies that a trial wave 
function forcing them to sit on the same site would not be an energy eigenstate. The same 
argument can be made the other way round. 



B. One-coloron wave functions 



In the following, we show by explicit calculation that momentum eigenstates constructed 
from |^/^) via (|2*5|) are exact eigenstates of i?su(3)- F° r simplicity, we choose a = b = y (an 
anti-blue or yellow coloron) and express j^^) through the corresponding wave function 



# 7 [zj; w k ] = JJ(r? 7 - Zj) #ota; «;*] = ^7*0, (27) 

i=l 

with given by © and M 1 = (N - 2)/3, M 2 = (N + l)/3. The momentum eigenstate 
h&^g) is then given by 

1 - 

w k ] = — J2(^) n * 7 [*i; ^1 ( 28 ) 

7=1 

(where n is shifted with respect to n in (|23jl by a constant depending on which momenta 
are occupied in the Slater determinant state). The momentum of (f^Hj) is 

4tt 2tt / l\ 
P=y-^U+3j, 0<n<Mi. (29) 

The energy is given by 



(30) 



e(p) 



allowed momenta 




2'tt 4'tt P 

T T 27r 



FIG. 3: One-coloron dispersion relation. 



where we have defined the one-coloron dispersion relation (see Fig. |3J) 



e(p) 



7T- 



7T 



(31) 



4 V 9 

For Mx < n < N, vanishes identically. 

To prove that *f> n represents an energy eigenstate, note that all terms of the Hamiltonian 
(jllj) yield the same expressions as for the wave function \&o except those terms containing 
e^e% and e^e T p. Hence, it is sufficient to investigate those further. 

To begin with, the polynomial ^ 7 defined in (12 7j) can be written as 
J2mLo(~ l) Ml ~ m rj™ S Ml (zi ■ ■ ■%!_„), where S M {z\ ■ ■ ■ z q ) is the sum over all possible 
ways to choose q coordinates out of z\, . . . , z M , e.g. S 3 (z 1 z 2 ) = z^z 2 + Ziz 3 + z 2 z 3 . Hence, 
we find in addition to (fT2j) - (fT5|) the two terms (note that deg z .\I/ 7 = deg Wfc \l/ 7 = N — 1) 



Mi / 1 

*«E 

i=i \ 



4 <9z 4 2 ^-f.z { - Zj dzi 



N-3 d \ , z; 



d 



dz. 



i=l fc=l 



2,- 



</V (32) 



Fourier transformation 

JV — 1" 
Mi I Mi 



of these yields 



ra n 



iV-1 



Ar Mi M% 

*-4*°EEE 



Furthermore 



E 



br rb 



-Vl> 



'/o 



A/i A/ 2 



7=1 i=l k- 



^ (zi - w k )(r]^ - Zi) 



W (33) 



i=l k=l 
N Mi M 2 



^EEE 



ZiWk^ 



Mi M2 

-p-r Zj - W k -p-r - Zj 

j~ % _ ^ 7M: w ' ~ w fc 



1 7 



(34) 



Mi 

n 



- w k 



The term on the right of (}3"4*)) is identical to a corresponding term in the ground state 
calculation, whereas (|H5|) taken together with the second term in (fHHJ) gives 



--n(n + l) + -M 1 (M 1 + l) 



where we have used 
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Theorem 3 For Mi, M 2 G N, Mi < M 2; z±, . . . , z Ml € C distinct, w±, . . . , % 2 G C distinct, 
and < n < M\ the following is valid: 



Mi M 2 

EE 

i=l k=l 



Zi - W k 



Mi M 2 



- Zi 



l^k 



Wi - w k 



N 



Mi 

7=1 j^t 



--n(n + l) + -M 1 (M 1 + l) 



/V 



Mi 



(36) 



7=1 



i=l 



The proof of Theorem is given in App. |Ej Collecting the terms obtained above, we find 
that \l/ n is an exact eigenstate of the Hamiltonian (0) with energy (|30|). 



C. Quantum numbers of colorons 

In the preceding section we have shown that the elementary coloron excitations are con- 
structed by annihilation of a particle of color a from an overall color singlet |^sb +1 ) before 
Gutzwiller projection. Since they are hole-like excitations, they transform according to the 
representation 3 conjugate to the fundamental representation 3 of the original particles on 
the sites of the chain. This can also be seen by acting with the total SU(3) spin generators 
(|21) on the wave function (|27|). Our result agrees with results obtained for the spectrum of 
the SU(3)i Wess-Zumino-Witten (WZW) model by Bouwknegt and Schoutens 3Q- 

It is straightforward to read off the quantum or exclusion statistics j^J \l\ of color- 
polarized colorons {i.e., colorons of the same color). Consider a chain with iV = 3M — 1 sites 
and a single yellow coloron. According to (JHJ), (J23J), and (|25j) . there are as many single particle 
orbitals available to the coloron as there are blue particles in the Slater determinant state, 
that is, M. If we now were to create three additional yellow colorons, the Slater determinant 
state would have to contain three more particles, one of each color. This implies there would 
be one additional orbital, while the three additional colorons would occupy three orbitals, 
meaning that the number of orbitals available for our original coloron would be reduced by 
two. The statistical parameter is hence given by g = 2/3. The fractional statistics manifests 
itself further in the exponents of the algebraic decay of the form factor of the dynamical 
structure factor [3, |^ as well as in the thermodynamics of the model j^, 2^]. A similar 



exclusion statistics exists in the conformal field theory spectrum of WZW models 



28, 



The exclusion statistics among colorons of different colors as well as the state counting for 
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SU(n) (n > 3) spin chains in general is highly non-trivial and will be subject of a future 
publication jsj]. 



V. TWO-COLORON STATES 

For N = 3M — 2, there are at least two colorons present. Two yellow colorons localized 
at lattice sites r/ 7 and rjs are described by the wave function 



(37) 



i=i 



where M\ = (N — 4)/3 and M2 = (N + 2)/3. Momentum eigenstates are once again 
constructed by Fourier transformation 



1 N 

tfmnfe; W k ] = Yj^TimY^-ysiZi] W k ], < 71 < m < M X . 



(38) 



7*5 



In analogy to the SU(2) HSM we neither expect these states for fixed m and n to be energy 
eigenstates, nor states with different sets of quantum numbers (m, n) with the same m + n 
(and hence the same total momentum) to be orthogonal to each other. Energy eigenstates, 
however, can be constructed as follows. 

As for the one-coloron states, application of the Hamiltonian generates two contributions 
in addition to those familiar from the the ground state calculation. First, for the terms 
containing e^ g e| b we obtain ([3*2*)) with ip 1 replaced by -?/V- The ^ rs ^ term in the result- 
ing expression can be treated in analogy to the two-spinon states in the SU(2) HSM [35]. 
Specifically, we find 

-Zi— ) </v = M l {2M l -N + 2)^ s 



2s I 2 Zi dzf + ^ Zi - z ,dz., 



i=l 



+ 



2 &_ 

77 w 



°7 

N-3 



'dm 



2 

d 



V>7<5 



<9r/ 



5^5 



where we have used (II 6|) three times. Applying this identity to the momentum eigenstates 
(jHHJ), we obtain 



-Mi (Mi + 2)\P TOn + 



~M X + 1 - m J m + ( -Mi + 1 - n J n — 



* 7 



^(m-n + 2£)^ m+£jn _ 
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where £ m = min(Mi — m,n), and we have used 

m — 7i. 

x + y 



(x m y n - x n y m ) 



1=0 



x-y 

The second term of (J3"2"j) applied to the states ^ mn yields 

1 



2j2 xm ~ £ y n+e ~ (x m y n + x n y m ), m>n. 



(39) 



N Ml M2 zKvTimT 



Zi - w k 



+ 

% - Zi T)5- Zi 



(40) 



7<5 i=l k=l 

The second contribution not familiar from the ground state calculation is given by terms 
arising in 



' N 

E 



br rb 



-\1> 



^0 \ ^ \ ^ \ ^ ZjW k 

+ /V2 _ 



A/i M 2 Mi M 2 

^ ^ ZjWfc -p-r Zj -W kT1 Wl-Zi 

N Mi Ah f ]m( , n Mx 

ZiWk{r] y ) [T)s) 



M 2 

Zj — W k "TT Wl-Zi 



TT fj ~ w fc TT 

^ 2 ^ tr (* - ^)(^ - *) }i - ^ U w < - 



(41) 



(42) 



JV Mj M 2 . , m( , n Ah Ah 

, W 0V^V^V^ ZiW k {f]^) [T] S ) T-T Zj - W k yj Wi - Zi 

+ m2^l^2^(r._ „„ - ^ 11 7T7 11 ^-^^ (43) 



■yS i 

N Ah Ah 



^0 v 1 ^ ZiW k (r] 1 ) m (r]s) n yr Zj - Wk jj Wi - Zj j 

N2 hhh <<h - *)(* - *> ii * - * i\ «* - «* 7i ' ( ' 

The first term (|41j) is identical to the ground state calculation. The two following terms (J42)) 
and (J43)) can be combined together with (J4Uj) and Theorem |3] to yield 



M 1 (M 1 + 1) - -m{m + 1) - -n(n + 1) 



(45) 



The last term (|4*4"|) can be simplified using 



Theorem 4 For M 1 ,M 2 G N, Mi < M 2 - I, z x , ■ ■ ■ , z Ml G C distinct, wi, . . . , ^Af 2 £ c 
distinct, and < n < m < Mi, and £ m = min(Mi — m, n) the following is valid: 

Mi M 2 A>h AI 2 N Mi 

E 5>» II II S IK* - - *> 



1 = 1 k = l JJ±l J lytk 



75 



N 



Mi 



75 



i=i 



5 - Zi 



Mi 



(46) 



■yS 



i=l 
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A proof of Theorem 0] is given in App. [FJ 

After collecting all terms, we obtain for the action of the Hamiltonian on the momentum 
eigenstates (J3S|) : 

vr 2 / 17 52 



\^mn) (47) 



#SU(3) |*mn> = "18 l^"^*^) l*™> 



3tt^ 
+ iV 2 



2 

(Mi — m)m + (Mi — n)n (m — n) 

3 



9 



y^(m - n + 2£) |# m+ ^_^) . 




This Sutherland-type equation [36, 37] shows first of all that when we act with the Hamilto- 
nian on the non-orthogonal basis state (J38|) with m = m , n = n , mo > no, we only obtain 
terms (J38|) with m > m and n < n with fixed m + n = m + no, but no terms with m < m 
or n > n . Consequently, states with m = Mi or n = (or both) are exact eigenstates. 
Furthermore, we can obtain all the other exact eigenstates at each fixed m + n (or fixed total 
momentum) by subsequently constructing an orthogonal basis of states starting from these 
exact eigenstates. (We have shown previously j3] that this method can be used to obtain 
all the spin polarized two spinon eigenstates of the SU(2) HSM.) 

With all the numerical constants in place, however, (J47j) provides much more information 
than the direction of the scattering. We can solve for the energy eigenstates directly through 
the Ansatz 

$mn = a T n ^m+i,n-i- (48) 

1=0 

Requiring 

Hsv(3) l^mn) = E mn |$ mn ) (49) 

then yields the recursion relation 



a 



5 )(m - n + £ - l)(m - n + 21) 

2 



mn V 3 ' ^ / V * „mn „mn 



£(m - n + £ + |)(m - n + 2£ - 2) 
i™ = —4/5) for the coefficients a™" as well as 



a<?_i, % = 1, (50) 



7T 2 / T 17 52 \ 3tt 2 
E mn = -— N - — + —) + 



2 



18 V N N 2 J N 2 
for the two-coloron energies. 

As we introduce single-coloron momenta according to 



(Mi — m)m + (Mi — n)n [m — n) 

3 



(51) 



4 2n ( \ 4 2it 

Pm = r~N [ m+1 h Pn =r~N l n+ i)' (52) 
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the energy ()51|) simplifies to 

4 7T 2 

E mn = E + -— + e(p m ) + e(p n ) (53) 

with Eq given by (fTU|) and the single-coloron dispersion relation e(p) given by ()31|) . The 
interpretation of (|52|) and (|53j) is as follows. First, the coloron excitations in the SU(3) 

) free, and the total energy is 



HSM are (like the spinon excitations in the SU(2) model 
simply given by the sum of the kinetic energies of the individual colorons. This conclusion 
is consistent with the asymptotic Bethe Ansatz calculation by Essler j^, who has shown 
that the coloron-coloron scattering matrix is trivial, i.e., a momentum independent phase. 
Second, the difference in the individual coloron momenta ()52|) is quantized as 

2tt (1 \ 
Pn-Pm = ^7\- + integer , (54) 



N \3 

i.e., the minimal momentum spacing between color-polarized colorons is 2/3 of the momen- 
tum spacing for fermions on a ring, 2n/N. We interpret this result as a manifestation of 
the fractional statistics of the colorons. This result appears to indicate that the spacing for 
particles with exclusion statistics g is given by 

Pn ~ Pm = (g + integer) . (55) 

This identification is obviously consistent with the familiar cases of bosons and fermions, 
and may even constitute the most direct manifestation of Haldane's fractional exclusion 
principle per se. 

The fractional momenta (|52|) may also be interpreted in terms of an effective change in 
the periodic boundary conditions (PBCs). The allowed values for the total momenta of the 
two-coloron states are those for conventional PBCs, 

Pn+Pm = j-jr ■ integer (56) 

where we have used the fact that iV + 2 in (J52j) is divisible by three. This is also required, 
as the wave function must come back to itself under a full translation around the ring. 
The allowed values for the difference in the momenta (J54j) . by contrast, are those of a ring 
threaded by a magnetic flux 27T-2/3. In a certain sense, we may interpret this flux as a phase 
acquired by the wave function as one coloron goes "through" the other. There are, however, 
several subtleties associated with this interpretation j^j. 
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VI. SU(3) SPIN CURRENTS 



In this section, we calculate the SU(3) spin currents, i.e., the eigenvalues of the diagonal 
rapidity operators A 3 and A 8 . To this end, we first rewrite them in terms of colorflip 
operators: 

TV 

(i ...((./ I l \ 

(57) 



A 



1 N 



Va+VP (> rb , £ bg gl 



-e rg e gr 



A 



y'3 //,, - //,; / !,„ ,!, 

A 2^ „ -„A e * 6 P 



4 ZA^-™ 



P "& P B» I p rg gr 



(58) 



The eigenvalue of A 8 is easier to evaluate. Using the Taylor expansion technique we obtain 

v^ 2 



m=0 



E 

,i=i 



^-(N-2)(M 1+ M 2 )*-^- 



' Mi 



m 2 



d * 



»7 



8=1 



5^ f-f (9w fc w fc 
fc=i 



(59) 



where the wave function has to satisfy deg^ . ^ < N — 1 and deg^ \I/ < N — 1 , and we have 
defined B m = J2a=i VaiVa + l)(Vot — l) m_1 . Evaluation of the latter yields B = N — 2, 
B x = -2, and S m = for 2 < m < N - 2 (see App. EJ. 

With (p)9"j) we find for the ground state, the one-coloron states ([28)1 . and the two-coloron 
states (jlHj) 



A 8 |^o) 

A 8 |^ n ) 

A 8 |$ mn ) 



0. 



N — 2 



V3( 
4 V 3 

(2N-8 



-2n |* n ), 



- 2m - 2n) |$ mn ) 



(60) 
(61) 

(62) 



4 V 3 

The eigenvalues of A 3 can be obtained either by explicit calculation or be determined as 
follows. The fundamental representations of the Yangian Y(sls) can be constructed 
from the representations of sl 3 (i.e., SU(3)) by the pull-back under the so-called evaluation 
homomorphism from K(sls) into the universal enveloping algebra £/(si3) of SI3. Explicitly, 
the ground state transforms according to the pull-back of the SU(3) singlet representation 
under K(sl3) transformations, the one-coloron states transform according to the pull-back 
of the representation 3, and the two-coloron states investigated in Sec. |V| according to the 
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FIG. 4: Weight diagrams of the SU(3) representations 3 and 6. 



pull-back of the representation 6 (as the two colorons are coupled symmetrically). For the 
fundamental representations of ^(skj) constructed in this way the ratio of the eigenvalues 
of A 3 and A 8 is equal to the ratio of the eigenvalues of J 3 and J 8 . Hence, as the states 
\if? n ) and |$ mn ) are SU(3) lowest weight states (i.e., they are annihilated by I = J 1 —iJ 2 , 
V~ = J A — iJ 5 , and U~ = J 6 — iJ 7 ; see Fig. HJ), the eigenvalues of A 3 are simply \f?> times 
the eigenvalues of A 8 given in (|HU|) - (|H2|) . 



VII. ASYMPTOTIC BETHE ANSATZ 



Let us now briefly compare our results with conclusions drawn from the asymptotic Bethe 
Ansatz (ABA). We wish to stress that the ABA is not an exact method. In the derivation 
of the ABA equations one assumes the existence of an asymptotic region in which the 
particles (or colorflips) are non-interacting ; an assumption which is not valid rigorously. 
In particular, the ABA equations have no direct connection to the exact eigenstates of the 
Hamiltonian [3] as it is the case for the standard Bethe Ansatz. It turns out, nevertheless, 
that the ABA reproduces the spectrum of the model, hence there is an indirect relation 
between solutions of the ABA equations and energy eigenstates. Let us now set up the 
formalism. The eigenstates of the SU(3) HSM are specified by two sets of pseudomomenta 
k\ and kf satisfying the ABA equations 3 

M 2 M 1 

Nk] = 27T/1 - tt sgn(^ - k]) + vr ]T sgn^ 1 - fcj), 
tt sgn(^ 2 - k)) = 2nl? + 7i ~ fcj), 

3=1 3=1 
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with two sets of mutually distinct integer (or half-integer) quantum numbers in the range 

<^(N + M 2 - M 1 -1), < ^(M 1 -M 2 - 1). (64) 

In ()63|) we have restricted ourselves to 1-strings, i.e., the fcj' 2 's are real numbers, which 
means that we can only describe color-polarized states (for the general case see 
The energy and momentum depend only on the first set of pseudomomenta, 

i=l 
M 1 

p = J2(kl + 7t) mod27r, (66) 

where k\ G [— ir, ir]. 

The ground state is obtained by arranging the I^'s densely around zero, i.e., by choosing 
M 1 = 2M and M 2 = M for a chain of length N = 3M. Excitations are created by 
introducing holes into these sets of quantum numbers j^, e.g. for N = 3M — 1 we obtain a 
single hole 1^ in the rank-2 quantum numbers by choosing M 1 = 2M — 1 and M 2 = M — 1. 
In all practical applications we have to take the thermodynamic limit, such that the ABA 
equations (|6*3|) turn into integral equations. We introduce pseudomomenta densities 

N+l K i N+i N 

as well as bare hole densities a^ 2 = Y^j 8(k — Aj' 2 ) which describe the excitations. Then the 
ABA equations (|53*jl turn into integral equations with (^-functions as integral kernels 

a\k) + al(k) = ^-a 1 (k) + a 2 (k), 

271 (68) 

a 2 (k) + al(k) = a\k)-a 2 (k). 

The latter are a special and distinguishing feature of the HSM. For the ground state we have 
a^j 2 = and (fF)HJ) is solved by <7q = N/3n and a 2 = N/Qn [39]. This shows the simple nature 
of the ground state, which is explicitly given by the Gutzwiller state ©• It is, however, not 

1 2 

possible to derive the Gutzwiller wave function from the ground state densities a ' as there 
is no direct connection between solutions of the ABA equations and the exact eigenstates of 
the Hamiltonian. 
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One-coloron states correspond to single holes in the rank-2 pseudomomenta, i.e., a\ = 
and a\ = 5(k — A 2 ), which implies 

a\k) = al- l -al a\k) = a 2 - jja 2 . (69) 

There is no "hole- dressing" in the surrounding pseudomomenta, or equivalently, the pseu- 
domomenta densities are only affected through trivial shifts as we introduce holes. This 
observation reflects once again the fact that colorons in the SU(3) HSM (as well as spinons 
in the HSM) are free j^. The densities yield the correct coloron momenta ()29j) and 
energies (|3(jp in the thermodynamic limit. Holes in the rank-2 pseudomomenta represent 
excitations which transform under 3 3^|, hence we find again that colorons possess comple- 
mentary colors We wish to stress, however, that this result cannot be derived using 
the ABA exclusively. First, it cannot be answered whether rank-1 or rank-2 holes represent 
elementary excitation, and second, the ABA is not an exact method in contrast to calcula- 
tions based on explicit wave functions. From (jHHJ) we can further read off that when adding 
one coloron, the number of available orbitals described by cr 2 (k) is reduced by 2/3. It is, 
however, not completely clear to us whether this can be interpreted directly as a derivation 
of the exclusion statistics of polarized colorons. 

VIII. GENERALIZATION TO SU(n) 

The results derived here generalize directly to the SU(n) HSM. Consider a chain with 
one particle per lattice site carrying an internal SU(n) quantum number which transforms 
according to the fundamental representation n of SU(n). The Hamiltonian is given by 

27r 2 N - e aT e T p 

fl su W = # EE|J_ p . ( 7 °) 

where the operator e a J annihilates a particle of flavor r at site i] a and creates a particle of 
flavor a at the same site. (Once again, r and o may be equal in ()70|) . In comparison to 
(J2J), we have omitted a constant.) The model again possesses a Yangian symmetry [ljj and 
hence is integrable. 

If we use a polarized state of particles of flavor n as reference state and label the coordi- 
nates of the particles of flavor cr, 1 < a < n — 1, by zf, 1 < i < M a , the wave functions jstj 
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n-1 M a n-1 M a M T n-1 M a 

*o[<i = n n« - 4) 2 n n n« - 4) n n < 

T =1 ? CT<T 1=1 j = l cr=l 1 = 1 

constitute exact eigenstates |2l| of the Hamiltonian (fTUJ) . For A" = nM, M CT = M, (1TT]) is 
the ground state of (fTUj) with energy 

7T 2 /n - 2 2n- 1\ 

£ »=-12(— A ' + — J' (?2) 

The total momentum of this ground state is p = (n — l)ixM mod 2tt. 
Localized SU(n) spinons are given by 

Mi 

^[<] = n^-4)^o[<] (73) 
1=1 

for N = nM - 1, Mi = M - 1, and M 2 = . . . = M n _ x = M. As hole-like excitations they 
transform according to the representation ii. Energy eigenstates are obtained by Fourier 
transformation of (|73|). 

1 - 

*»[*f ] = T7 < m < Ml (74) 



N , 

7=1 



Their momenta are given by 



n — 1 27r/ n — 1\ , . 

p = vrN - — m H mod 2tt, (75) 

F n A/\ 2n / ' v ; 

which fill the interval [— -, -] for n even and M odd, or the interval [— -, -] otherwise (either 
n odd or M even or both). The energy eigenvalues of the one SU(n) spinon states (|74JI are 
given by 

E n (p) = E% + + 6«(P) (76) 

with (see Fig. EJ) 



e n (p) = 4 



^2 



, p 2 I , if n even and M odd, 

4 \rr / 

— ( — — (p — 7r) 2 ) , otherwise. 
„ 4 \rr / 



(77) 



The statistical parameter for polarized SU(n) spinons is g — (n — l)/n. In the framework 
of the ABA, the spinons are represented by holes in rank-(n — 1) pseudo momenta. 
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a n even 



e(p) 



b) n odd 




I I 
2vr 




2tt 



FIG. 5: SU(n) spinon dispersion relations, a) n even. The allowed momenta fill the interval 
\tt — —, tt + -1 for M even and f— — , — 1 for M odd. b) n odd. The allowed momenta fill the interval 

Finally, let us add some remarks on the large-n limit of the SU(n) HSM. First, there are 
no simplifications in the calculations which yield the results presented in this section. When 
n grows, we obtain solely terms similar to the ones treated in the appendix. Second, the 
interval of allowed spinon momenta shrinks to zero. This does not mean that the number of 
spinon orbitals vanishes, as in the limit n — » oo the number of lattice sites N = nM — 1 has 
to grow and hence the momentum spacing tends to zero as well. Third, we find g — > 1 as 
n — > oo, meaning that the exclusion statistics between polarized spinons becomes fermionic 
in this limit. This does not imply, however, that spinons in the SU(n) HSM behave like free 
fermions in the large-n limit. Their momenta fill only a narrow interval and the individual 
spinon momenta are not simply given by integer multiples of 27r/N 



IX. CONCLUSION 



In conclusion, we have shown by explicit calculation that the elementary excitations of 
SU(n.), but in particular SU(3), HSM transform under the representation conjugate to the 
representation of the SU(n) spins on the chain. We have also shown that they obey fractional 
statistics, and only exist in one n th of the Brillouin zone. We strongly believe these features 
to be generic features of "antiferromagnetic" SU(n) spin chains. We have further endeavored 
to construct the polarized two-coloron eigenstates of the HSM explicitly, and have thereby 
shown that they are free and that the spacings in the difference of the individual coloron 
momenta p n and p m with p n > p m are given by p n — p m = (| + integer). We have 
interpreted this peculiar shift as a manifestation of the fractional statistics. This indicates 
that fractional statistics in one dimension manifests itself not only through a fractional 
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exclusion principle 24], but that there are Berry's phases associated with anyons 23j in 
one dimensions as well 40] . Furthermore, we look forward to an experimental realization of 

fin 

SU(3) spin chains in systems of ultracold atoms [12], [13J, and hence an observation of the 
complementary colors of the excitations. 
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APPENDIX A: GELL-MANN MATRICES 

The Gell-Mann matrices are given by 





f° 


1 






f° 


—i 






(l 
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A x = 


1 










i 













-1 






1° 





V 




\° 





V 




1° 





/ 




f° 





A 




f° 





A 




f° 


0^ 




A 4 = 











, A 5 = 











. A 6 = 





1 


1 




V 





V 




V 





J 




1° 


10J 





^ ^ 



A 7 



0-i 
i 



1 

71 



1 
1 
0-2 



\ 



/ 



They are normalized as tr (A a A fe ) = 25 a b and satisfy the commutation relations [A a , A b ] = 
2f abc X c . The structure constants f abc are totally antisymmetric and obey Jacobi's identity 



ya6c jcde _|_ jbdc jcae _|_ jdac jcbe 



0. 



Explicitly, the non-vanishing structure constants are given by f 123 = i, f 1Ar = f 2m = f 257 = 
j345 = _yi56 = _j367 = f y 2> / 458 = f 78 = iy/3/2, and 45 others obtained by permutations 
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of the indices. The SU(3) spin operators are expressed in terms of the colorflip operators by 



2^ a p 6 

a=l err 

APPENDIX B: EQUIVALENCE OF © AND © 

Let us begin by exploring a general consequence of particle-hole symmetry. Consider a 
chain of N sites. The wave function of L fermions with momenta in the interval I\ is given 
by 

$!K ...,u L } = (0| c ui ■ ■ ■ c UL n c\ |0>. (Bl) 

qeh 

This wave function can be expressed in terms of N — L fermions whose momenta fill the 
interval I2 = [— 7r, 7r]\Ji, 



U L ] = (<0| II C„ ^ - - -C+ x |0> 



= sign • ( (0| Y[ c q c Vl --- c VN _ L Y[ c\ Y[ c\ |0) j 

^ qeh qeh qeh ' 

= sign-m Vl ,...,v N - L }, (B2) 

where the coordinates v\, . . . ,vn-l are those coordinates on the chain which are not con- 
tained in the set U\, . . . , ul, 

$ 2 [vi,...,v N - L ] = (0\c Vl ---c VN _ L IJ c+ 10), (B3) 

qeh 

and the overall sign depends on the number and positions of the fermions. 

(jB2|) enables us to show that the ground state wave function is given by with M\ = 
M 2 = M. For N = 3M, the wave function of the Slater determinant state |^sd) is ( see 
Fig. Hi) 

^ST>[Zi]W k ;u m ] = $i[z,]$iK]fi[ii m ], (B4) 
where the interval I\ is chosen to be [— 7r/3, tt/3] and hence 



M „ , M 

M-l 



i=l i<j 

The prefactor in ()B5|) shifts all momenta from the interval [0, 27r/3] to I\. The ground state 
of the SU(3) HSM is obtained by Gutzwiller projection, |\& ) = Pq I^sd), which enforces 
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that all coordinates Zi, Wk, and u m are distinct. Hence, we can express the coordinates u Ti 
in terms of the z^s and w^s using ()B2|) . 



$i[u m ] = sign ■ ^* 2 [zi,w k ], 

where the wave function of the filled band with momenta in the interval I2 = [— tv, is 
(up to a sign) given by 

M M M M M M 

M+l A/+1 

* 2 h,w fc ] = J_J_2;. 2 J_J_w fc 2 J J_(2i — U(wfc — wi) 11 — ( B6 ) 

i=l fc=l i<j fc</ 6=1 fc=l 

The ground state wave function is hence given by 

M M M M MM 

i=l k=l i<j k<l 1=1 k=l 

where we have used the fact that the coordinates are located on the unit circle, e.g. z* = 
Using 



(7 7 \ 2 M 1 M 1 

12 _ %J TT— — TT 

11 r r 11 jw-r 



* Z 3 1 'II 11 M— J ' 

. , . ZjZj . Z- 
i<j J i=l 1 



as well as 

we obtain from ()B7|) 

M MM M M M M M 

%[zi, w k ] = n Z f M n W f M - - n - ^) n * n 

i=l fc=l i<j fe<i i=l fe=l i=l k=l 

As the coordinates z$ and w k are located at the lattice sites i] a = exp (^pa) and iV = 3M, 
the prefactor reduces to one and \E r o[zi,u'fc] takes the form (JSJ). 



APPENDIX C: USEFUL FORMULAS 

Some of the results presented in this appendix can be found 
I. 

N N 



in 0. 



^=1, £c=^<w, n^=(- i ) 7v ~ i - ( ci ) 



o=l a=l 
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FIG. 6: Contours for integrations 

2. 



N 



]J(z- Va ) = z N -l. (C2) 



a=l 



3. 



Proof: Rewrite (jC3|) as 



N N 



with / a polynomial of degree N — 1, /(0) = AT, and /(l/^ 7 ) = n^Ji^l — Vp))! = 
1, . . . ,N. Since the polynomial / — /(l) has de gree N — 1 and iV different zeros, it 
vanishes by the fundamental theorem of algebra [42]. 



4. The previous statement implies 



N 



7V-1 



5. 

JV-1 



-ri m AT _|_ 1 

V — ^- = — — -m, l<m<N. (C5) 

^ ?]rv - 1 2 



Proof: Using Cauchy's theorem 42] for the function 



z JV (2; — 1) 



with the contours drawn in Fig. El yields 



N-l m , iV-1 



. 7U - 1 27T2 ^ 



^?7a-l 27Ti Z*(z- 1)3-77, 
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N f z m - x 

dz 



2m Jc (z - l)(z N - 1) 

n r z m - x 



2m J c , (z-l){z N - 1) 
N+ 1 



dz 



— m, 



2 

where we have used the theorem of residues 12] in the last step. 



N-l 



Proof: We have 

N-l _, JV-1 



m+1 



In - 112 w , 



^ ha - IP ^ fa, - I) 2 

N 



2m Jc (z - l) 2 (z N - 1) 
NT z m 



dz 

2m J c , (z-l) 2 (z N -I) 

N 2 - 1 m(JV-l) m(m-l) 



12 



+ 



7. For 



2V-1 



a=l 



we have: A = (JV - 1)(JV- 5)/12 by (Jd, A x = -(iV-3)/2 bv (jHH|) . and A 2 = 1 
by Va = N5 m o. Furthermore, 

7V_1 m ~ 2 An 2\ 

a=l fc=0 ^ ' 
m-2 / \ AT 

= E( m ; 2 )(-ir- l (i-E"r) 

fc=0 ^ ' a=l 

= (f) (~ l ) n ~ k = °> 2 < m < iV — 1. 

fc=0 ^ ' 

as the sums of the binomial coefficients of even sites and odd sites equal each other. 



8. For 

N-l 



Of=l 



we have: B = N - 2, B t = -2, B m = for 2 < m < N - 2, and 5tv-i = iV. 
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APPENDIX D: PROOF OF THEOREM Q] 



As (J22J) is symmetric in the Zj's, we prove the theorem by showing that the LHS of 
is independent of Z\ and vanishes identically. To this end, consider the auxiliary function 

*(*i) = £ #- !) — + E — -^4 • ( D1 ) 

Y[ j=2 {Zj - Zi) ^ (Zi - Zi ) \\ j+i {Zj - Zi ) 

for fixed and distinct z 2 , ■ ■ ■ , Zm £ C. 

First, we show that ()D1|) is entire (i.e., analytic on C). We can restrict ourselves to 
showing that is analytic at Z\ = z 2 . (jDl|) can be rewritten as 



®{Zl) = jj h <P(2l, 



where $ is analytic at Z\ = z 2 . Expanding the first term yields expressions like 

1 



-(^ 2 m -«), m,nGN, 

Z\ — z 2 

which can be shown to be analytic at Zi = z 2 by using 

Second, as the degree of the denominator in (jDl|) is strictly greater than the degree of the 
numerator, vanishes for \z\\ — > 00. Therefore, by Liouville's theorem j3] $ vanishes 

for all Z\ 6 C, which proves the theorem. 



APPENDIX E: PROOF OF THEOREM |3] 



We prove the theorem in three steps. First, we show that the LHS of (fMJ) does not 
depend on the WkS. Second, we eliminate the w^s. Third, we verify the remaining equation 
by explicit calculation. 

For the first step, note that (|36J) is symmetric under permutations of the Wk's. It is 
hence sufficient to show that (JHljj) does not depend on w\. Consider for fixed and distinct 
zi, . . . , zm 1 ,w 2 , . . . , wm 2 £ C the auxiliary function 



h - w x 

M 2 



Mi M 2 

n Z 3 ~ W l TT W l ~ Z l 



Zj — Z\ 



1=2 



Wl - Wi 



M 2 



k=2 



Mi M 2 

nzj - w k yj wi - Zi wi- z\ 
— z l 

3 _ 2 Z 3 ~ Z l t/k Wl ~ Wk Wl ~ Wk 



(El) 
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We will show that $ is an entire function and bounded for \w±\ — > oo. Then Liouville's 
theorem |42| implies that <l> is a constant, i.e., independent of W\. 

To see that $ is entire note that (|Klj) is obviously analytic for W\ E C\{z 1 , w 2 , ■ ■ ■ , %J. 
To investigate the point W\ — z±, we rewrite the first line of (jKl|) as 

M 2 



Zi - w x 



Ml 



1 + 



Z 3 



n 



Zl - Wl 
Wi - Wi 



Zl 



j=2 ^ ■> x ' 1=2 

which shows that $ can be analytically continued to w± = z\. Furthermore, as (|E1J) is 
symmetric in w 2 , ■ ■ ■ ,wm 2 , we may restrict ourselves now to showing that $ is analytic at 
W\ — w 2 . For w 2 = this is obviously the case. For w 2 7^ we write $ as 



$00 



Wi 



Mi Ah 
TT Zj - Wl yr Wi — Zj _ W 2 - Z 1 

Zi - Wi jj^ Zj - Z 1 Al Wj - U7l W2 - Wl 



+ 



w 2 



Ml 



II II 



M 2 



^1 - ^2 " ^ 



Wi - Zi Wl - Zl 



where <£> is analytic at W\ = w 2 . Now, $ — $ can be rewritten as 

1 Mi 1 M 2 



1 



: 'l f = 3 ( w f _ Wi)(w t - W 2 ) W X - W 2 



(Zi - Wi){Zi - W 2 ) " z. 

Mi M 2 Mi M 2 

wi(zi - w 2 ) 2 Y\[zj - wi) Y[(wi - w 2 ) - w 2 (zi - Wi) 2 \\{Zj - W 2 ) Y\i w i ~ w l) 



1=3 



3=2 



Expanding the square brackets in (jK2|) leads to terms like 

1 



J=2 



1=3 



(E2) 



(w?w?-w?w%), m,neN, 
Wi - w 2 

which can be shown to be analytic at Wi = w 2 by using J3UJ). This concludes the proof 
that $ is an entire function of w\. Furthermore, $ tends to a constant for \wi\ — > 00 as 
in the first term of (jEljl . the degree of the denominator is greater or equal than the degree 
of the numerator, while in the second term the Wi-dependence cancels out for \wi\ — > 00. 
Therefore, $ is constant by Liouville's theorem j^]. As we can replace Zi by any Z{ in ()E1|) . 
the LHS of (HSU does not depend on the w^s. 

In the second step, we eliminate the w^s from (J36|) . In order to do so we choose Wj = Zj 
for 1 < j < Mi, and w k = Re itpk with distinct values ^ k G [0, 2vr] for Mi + 1 < k < M 2 . We 
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now have to study three different terms in (JHUj) . First, the terms with i = k yield as R — > oo 

Mi N Mi 

-E^E^rn^ - ^)- 

i=l 7=1 M« 

Second, the terms with i ^ k,k < Mi yield in the same limit 



Mi Mi 2 N M i 

1=1 k=tt 7=1 j^i 



7 ^iJ- 



Third, the terms with Mi < k vanish in this limit. Hence ()36|) reduces to 

Mi N Mi Mi Mi 9 N M i 



Mi Mi 2 N 

E * E(^) n - **) + E E 737 E(^) n Ilfo - < ) 

2 = 1 7=1 i — >i-Ai ^ 3 



/7 i / / / J I I V77 

kjti i=l jjti 1 7=1 k=ii 

N Mi 

7=1 8=1 



(E3) 



with /(n) = -n(n + l)/2 + Mi (Mi + l)/2. The second term in the LHS of (JE3J) can be 
written as 



AT Mi 9 

zf 



M-i 



E^)" E 7T7^ - *i) II - z *) 

7=1 Mi < j 



*M*,i 



N 

E 

7=1 



Mi Mi 



Z%Zj 



Mi 



Mi 



Mi 

n - **) 



(E4) 



where we have used 



Mi ^ 2 j M i 

E TTT^ - *i) = o E + ^) ~ ^'1 • 

Zi Zr 



Mi ^ J Mi 
In the third step, we complete the proof by explicit calculation. For this we use 



N Mi 
7=1 i=l 



l) M ^S M ^ Zl ...z Ml . 



(E5) 



where S Ml (zi • • • ZMi-n) is the sum over all possible ways to choose Mi — n coordinates out 
of Zi, . . . , zmi- Then, using (jE4|) we write (jE3j) as 



Mi 



yZ\ ■ ■ ■ ZMi-n-l, 



E^ Ml_1 

i=l 

1 Mi 

+ 77 ^i^j 2 (^1 - - " ^Mi-n-2, 



Mi 



Mi 



Mi 



no Zi,Zj 



Z\ ■ • • ZMi-n-l, 



f{n)S M ^z 1 ---z M ^ n ). 



no z-i,Zj 



(E6) 
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With the relations 

Mi 



i=l 
Mi 



Zi S Ml 2 {z\ ■ ■ ■ ZM\-n-\. 

Mi 



no Zi,Zj 



(M 1 -n)S M ^z 1 ---z Ml -n), 
{M 1 -n)nS Ah {z 1 ---z Ml ^ n ), 

(M x - n)(Af! - n - 1) S Ah ( Zl ■ ■ ■ z Ml -n), 



we can verify (|E6|) and deduce f{n) = —n{n + l)/2 + Mi (Mi + l)/2. This concludes the 
proof of the theorem. 



APPENDIX F: PROOF OF THEOREM H 



Consider for fixed and distinct z\, . . . , zmi,w 2 , ■ ■ ■ , % 2 G C the auxiliary function 

Mi Mi 

j=2 ■? 1 1=2 

M 2 Mi M 2 

+ Tw k T] Z] ~ Wk T] Wl ~ Zl ■ Wl ~ Zl - 
f 2 7=2 z i~ Zl U Wl ~ Wk Wl ~ Wk 



(Fl) 



In analogy to the proof in App. |E| it can be shown by using ()39|) that (jFlj) is entire, and 
that $ tends to a constant for \w%\ — > oo. By Liouville's theorem 42] $ is constant, and 
hence is independent of the w^'s. We can choose the w k s as Wj = Zj for 1 < j < Mi, 
and w k = Re iipk with distinct values p k G [0, 2vr] for M x + 1 < k < M 2 . Then, the LHS of 
(EEoT) simplifies to 



Mi AT Mi 

E z " E^r^r II(% - - *,■)■ 

Using (jE5j) we can rewrite (|4*Uj) as 



(F2) 



Mi 



Z\ • • • ZMj-m-1. 



i=l 



S hh -\ Zl ---z Ml - 



•n-l. 



no z., 



(Mi - m) 5 Ml ZMi-m) S Ml («!••• ^Mi-n) 

£ 11 1 

J^( m _ „ + 2 £) 5 Ml (Zi • ■ ■ ZMi-m-^) $ Ml (Zf- Z Ml . 



(F3) 
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As all terms in (jF3|) are completely symmetric in the z^s, we simplify the notation by 
imagining to have written an operator S R!l around each term and then order all terms inside, 
i.e., we replace 

S M {z\ ■ ■ ■ z m ) — > ( 1 Z\ ■ • ■ z v 



m / 

The remainder of the proof makes use of the auxiliary theorem 
Lemma 5 

k. 



S"(z 1 -z m )S«(z 1 -z n )^ g (Z)(T)( M n-Z) lk ' m + n - k] ' (F4) 



where k min = max(0,m + n — M), k max = min(m, n), and [k, m + n — k] 



2 2 
%\ ' ' ' %k ^k+1 ' ' ' Zm+n—k ■ 



Proof: Replace first the LHS as 



S M (z\ ■ ■ ■ z m )S M (z\ ■ ■ ■ z n ) — > ( ) z\ • ■ ■ z m S M (zi ■ ■ ■ z r . 



In order to obtain the contribution z\ - ■ ■ z\, k coordinates out of Z\, . . . , z n have to match 
k coordinates out of z\, . . . , z m , which yields (T) terms for /c min < k < k max . The remaining 
n — k coordinates out of z±, . . . , z n match z m+ i, . . . , z m+n _ k , which yields ( A ^I™) terms. □ 

Now, by application of this lemma the first term on the RHS of (|F3|) simplifies to 
. A ^ m /MA /Mi - m\ ( m \ ri 

k — k m 

where k m = max(0, Mi — m — n) and we have used n < m. The second term on the RHS of 
(lF3l) reads 

v^V^ „n f Mi \ (M 1 -m- A / m + £ \ rj 

£=1 fc=fc m \ 1 / \ / \ l / 

In this equation we change the order of the summations according to 

^ m Mi-m-i Afi-m-1 Mi—m-k 

E E - E E ■ 

^ — 1 — /c 2n ^ — ^m ^ — 1 

On the LHS of (jF3|) at least is present and hence we find 

/Mi-A /Mi-m-l\ / m \ r , „ 

31 



where k' m = max(l, Mi — m — n). By considering each k separately, ()F3|) simplifies to 

k M\ — m 

(Mi - m - k)\(M 1 -n-k)\ ~ [M 1 -m- k)\(M 1 - n - k)\ 

Mi-m-fc „ (F5) 

(Mi - m-fc-^)!(Mi - n- + 

For k = Mi — m, (|F5|) is trivially satisfied. For max(0, Mi — m — n) < k < M\ — m — 1, 
()F5|) reduces to 

m n — m + 2£ 



which is easily verified. This completes the proof of Theorem |3J 
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